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We investigate the properties of mixed stars formed by hadronic and quark matter in /^-equilibrium 
described by appropriate equations of state (EOS) in the framework of relativistic mean-field theory. 
We use the non-linear Walecka model for the hadron matter and the MIT Bag and the Nambu- 
Jona-Lasinio models for the quark matter. The phase transition to a deconfined quark phase is 
investigated. In particular, we study the dependence of the onset of a mixed phase and a pure 
quark phase on the hyperon couplings, quark model and properties of the hadronic model. We 
calculate the strangeness fraction with baryonic density for the different EOS. With the NJL model 
the strangeness content in the mixed phase decreases. The calculations were performed for T = 
fT^ , and for finite temperatures in order to describe neutron and proto-neutron stars. The star properties 

■ are discussed. Both the Bag model and the NJL model predict a mixed phase in the interior of the 

' star. Maximum allowed masses for proto-neutron stars are larger for the NJL model (~ 1.9 Mq) 

CN I than for the Bag model (~ 1.6 Mq). 



< 



X 



PACS number(s): 2l.65.+i, 21.30.-x, 95.30.Tg 



I. INTRODUCTION 

Landau predicted the possible existence of a neutron star after the neutrons were discovered by Chadwick in 1932. 
, In 1934, it was suggested that neutron stars were formed after a supernova explosion, which happens when the core 
' of a very massive star undergoes gravitational collapse. The first supernova explosion was registered in 1054 by 
, the Chinese. The Crab Nebula in the Taurus constellation is the remnant of this explosion. Recently, a supernova 
' explosion was observed in the Magellanic Cloud 170,000 light years from the Earth. Once the gravitational collapse 
[ of a massive star with mass of the order or larger than 8 solar masses takes place, a proto-neutron star can be 
-4-^ . formed. Several different stages may happen during the evolution process The proto-neutron stars are known 

,_!h ' as evolutionary endpoints and they slowly cool down to form a neutron star, which is a stable and cold compact star. 
^ , The structure of compact stars is characterized by its mass and radius, which in turn are obtained from appropriate 



equations of state (EOS) at densities about one order of magnitude higher than those observed in ordinary nuclei. At 
these densities, relativistic effects are certainly important. 

In this work we investigate the equation of state (EOS) of warm, /3-equilibrium hadron/quark matter and apply it 
to determine the properties of mixed stars consisting of hadron matter with hyperons and quark matter. In particular 
we investigate the mixed phase with hadron and quark matter and search for the possibility of existence of a pure 
' quark matter core inside compact stars. The calculations are performed for finite and zero temperature in order 
to describe proto-neutron and neutron stars, respectively. We consider only the phase after deleptonization, when 
neutrinos have already diffused out. 

For hadron matter the relativistic non linear Walecka model 0,13 with the inclusion of the baryonic octet is used. 
For sufficiently high densities it is expected that the formation of hyperons in a neutron star is energetically favored. 
Experimental constraints obtained from hyper-nuclei give estimates of the hyperon-nucleon and hyperon-hyperon 
interactions, which impose restrictions on the expected properties of neutron stars @ . The appearance of the strange- 
baryons softens the EOS and lowers the maximum mass of a stable neutron star Another possible source of 
softening of the EOS is the onset of quark matter. According to the hyperon couplings and the properties of the 
quark model, this may occur for densities lower or higher than the one corresponding to the onset of hyperons 1]. For 
quark matter, the MIT Bag model ^] and the Nambu-Jona-Lasinio model f7| with three quark flavors were chosen. 
The NJL model contains some of the basic symmetries of QCD, namely chiral symmetry. It has been very successful 
in describing the low lying mesons and predicts at sufficiently high densities/temperatures a phase transition to a 
chiral symmetric state [1,13, 0, 01 ■ However, it is just an effective theory that does not take into account quark 
confinement. This rises no problems because we only use the NJL model to describe the deconfined phase. At 
sufficiently large baryonic densities the quark structure of the hadrons gives rise to a phase transition into quark 
matter. This conversion takes place at densities a few times the nuclear matter density 0. The Gibb's criteria are 
enforced in obtaining the coexistence phase. We also check whether the existence of the phase transition depends on 
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the choice of parameters for both hadron and quark matter equations of state, as pointed out in [l^]. The complete 
EOS is buih from the hadronic EOS at smaU densities, the mixed EOS for intermediate densities and a quark matter 
EOS for higher densities. 

One of the aims of the present work is to understand the importance of chiral symmetry for the description of the 
quark matter. There have been contradictory resuhs presented using the NJL model 13. 14iJ. namely with respect to 
the possibility of existing a mixed phase of hadrons and quarks inside a neutron star. We also study the role of the 
hyperon coupling constants in the appearance of the phase transition to quark matter and the change of strangeness 
composition with density. Mixed star properties are obtained by solving the appropriate equations for temperatures 
up to 30 MeV. 

The present paper is organized as follows: in section II both models used for the quark matter are described 
and in section III the non linear Walecka model with hyperons is reviewed. In both sections the conditions for 
beta equilibrium and charge neutrality are discussed. In section IV the mixed phase is implemented, the results are 
shown and discussed and in section V the properties of compact stars are computed. Finally, in the last section, the 
conclusions are drawn. 



II. QUARK MATTER EQUATION OF STATE 
A. The Nambu-Jona-Lasinio Model 



We describe the quark phase within a model with chiral symmetry, the SU(3) NJL model which includes scalar- 
pseudoscalar and the 't Hooft six fermion interaction that models the axial U{1)a symmetry breaking. The NJL 
model H, 0, 0, 0| is defined by the Lagrangian density 

8 

L = qiiYd^ - m)q+gs Y,[{qV qf + [qil^X'^qf] 

a=0 

+ to {det[gi(l + 75)gj]+det[9, (1-75)9^]}, (1) 

where q = (m, d, s) are the quark fields and Aq (0 < a < 8) are the U(3) flavor matrices. The model parameters 
are: m — diag (m„ , , ), the current quark mass matrix {md = the coupling constants gs and gu and the 
cutoff in three-momentum space, A. 

The set of parameters is chosen in order to fit the values in vacuum for the pion mass, the pion decay constant, the 
kaon mass and the quark condensates. We consider the set of parameters [ll|,[l^: A — 631.4 MeV, gs = 1-824, 
go^^ = —9.4, rriu = md = 5.6 MeV and nis — 135.6 MeV which where fitted to the following properties: m.^ — 139 
MeV, U = 93.0 MeV, mx = 495.7 MeV, !k = 98.9 MeV, (mm) = (Jd) = -(246.7 MeV)^ and (ss) = -(266.9 MeV)^. 
The frequently used set A = 602.3 MeV, = 1-835, go = -12.36, to„ = = 5.5 MeV and m, = 140.7 MeV 

|16| gives similar results but since it has a smaller cut-off and we want to study quark matter at high densities we 
have preferred the first set. 

The thermodynamical potential density is given byfi = £ — T S — ^^iNi — fig, where the energy density is 

£ = -2A^c 2^ / {n^- -n,+ )e{k -p ) 

- 2gs -'2 9D{uu){dd){ss) -£o (2) 

i—u,d,s 

and the entropy density is 

S = -2N, J2 y7^^(A2-p2){[n,+ln(n,+ ) + (l-n,+ )ln(l - + (3) 

i—u,d.s 

In the above expressions Nc — 3, T is the temperature, fii (Ni) is the chemical potential (number) of particles 
of type i, and £0 and are included in order to ensure £ = ^ = in the vacuum. This requirement fixes the 
density independent part of the EOS. The ground-state of the system is described by the density matrix given 
by / = diag (/„ , fd , fs) with 

/, = -[/(n,_-fn,+ ) - ^-^^i P(n._-n,+ )]0(A2-p2), (4) 
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where / is the identity matrix, rS-^'^ are the Fermi distribution functions of the negative (positive) energy states, 
■"'^^^ = [1 + exp(=F(/3 ± /ii)))]"-'^, i — u,d,s. In the last equation (3 = l/T, Mi is the constituent quark mass. 



The quark condensates and the quark densities are defined, for each one of the flavors i — u,d, s, respectively as: 

{q,q,)^~2N, j - n,+ ) e{K^ ~ p'), (5) 



= (g,t q^) = 2N, J -0^ + n,+ - 1) 9iA^ - p^). (6) 

Minimizing the thermodynamical potential Q with respect to the constituent quark masses Mi leads to three gap 
equations for the masses Mi 

Mi = nii -igsiqiQi) - 2gD {Qj qj){qkqk) , (7) 

and cyclic permutations of i, j, k. 

We introduce an effective dynamical Bag pressure 0, 0| , 

B = 2N, ^ J ^(^^^TmI-^^T^'^ 9{A^-p') 

i—u,d.s 

- "^gs ^ {qiqi)'^ - igo {uu){dd){ss). 

i—u.d,s 

In terms of this quantity the energy density (O takes the form 

£ = 2N, J2 j 0^Vj^^TM,{n,+ ~n,^ + l)6{K^ -p") +B,ff, B,ff = Bo-B, (8) 

i—u,d.s 

where Bq — -Bp^=p^=p,=o- Writing the energy density in terms of i?e// allows us to identify this contribution as a 
Bag pressure and establish a relation with the MIT Bag model 0,^3 discussed in the next section. 

We point out that the NJL model is only valid to describe the quark phase as far as the momenta of the quarks are 
smaller than the momentum cut-off A. 

B. The MIT Bag Model 

Quark matter has been extensively described by the MIT Bag model Q. In its simplest form, the quarks are 
considered to be free inside a Bag and the thermodynamic properties are derived from the Fermi gas model. The 
energy density, the pressure and the quark q density are respectively given by: 

£ = 3x2^ / 7^\/p' + (A+ + /?-)+ Bag, (9) 

q—u,d,s 

P=\T.f dp (A+ + /,-) - Bag, (10) 

P. = 3 X 2! ^(/,+ -/,„), (11) 

where 3 stands for the number of colors, 2 for the spin degeneracy, ruq for the quark masses. Bag represents the bag 
pressure and the distribution functions for the quarks and anti-quarks are the Fermi distributions 

/«± = l/(l + exp[(eTA^g)/T]) , (12) 
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with iiq being the chemical potential for quarks and anti-quarks of type q and e = yp-^ + yn^. These equations 

were obtained for finite temperatures. For T = 0, the expressions can be read off the above ones by eliminating the 
antiparticles and substituting the particle distribution functions by the usual step functions. 

We have used m„ = = 5.5 MeV, rus = f 50.0 MeV and Bag = (180 MeV)"* or Bag = (190 MeV)''. If m„, rud 
and TOs are chosen as in the NJL model, the behavior of the properties of interest are not altered, since they are more 
dependent on the Bag pressure than on small differences in the quark masses. 



C. Quark matter in beta equilibrium 

In a star with quark matter we must impose both beta equilibrium and charge neutrality j^l]. We consider the 
stage after deleptonization when entropy is maximum and neutrinos diffuse out. In this case the neutrino chemical 
potential is zero. For /^-equilibrium matter we must add the contribution of the leptons as free Fermi gases (electrons 
and muons) to the energy and pressure. The relations between the chemical potentials of the different particles are 
given by 

Ms = = M« + Me, Me = Mai- (13) 

For charge neutrality we must impose 

Me + Mm = ^(2mu - Pd- Ps)- 
For the electron and muon densities we have 

Pi = 2 j ^{fi+~fi^), l = e,p (14) 

where the distribution functions for the leptons are given in ea. H12|) by substituting q by I, with pi as the chemical 
potential for leptons of type I. At T = 0, equation H14|l becomes simply pi = k^i/iir^ . 

In figurefflwe plot the effective bag pressure B^ff for several temperatures and for quark matter in beta equilibrium. 
As discussed in [iJl for T = MeV around 161-163 MeV there is a plateau between p = Spo and 5mo- This corresponds 
to partial chiral symmetry restoration for quarks u and d but the chemical potential for quarks s, Ms, is still lower 
than rris. For higher densities, Ms > '^s and the mass of the quark s goes slowly to its current quark mass value as 
density increases. For high densities, if the current quark masses were zero, the bag pressure i? — > and the quarks 
would behave as a gas of massless non- interacting particles inside a large MIT Bag with a bag constant Bq [13 • At 
finite temperature, the value 161 MeV is only reached at higher densities and the plateau slowly disappears. This 
behavior is due to the fact that the mass of the quarks decrease slowlier with p and the quark s exists for p^ < In 
fig|21we show for different temperatures the fractions of m, d and s quarks, Yi = pi/{Nc p), as a function of p, the total 
baryon density. It is clear the effect of temperature on the appearance of the quark s. We will see in section IV that 
the behavior of the bag pressure with density and temperature determines the onset of the mixed phase corresponding 
to the coexistence of hadron and quark matter. 



III. HADRONIC MATTER EQUATION OF STATE 

An extension of the NLWM jj] is the inclusion of the whole baryonic octet (n, p, A, S+, YP , S^, S^, S°) in the 
place of the nuclconic sector. The presence of baryons heavier than the nucleons is expected in matter found in the 
core of neutron stars. The spin |A and the Jl" hyperon should appear at even higher densities and hence are not 
included. The inclusion of other meson fields describing the /o(975) and (^(1020), which are a scalar and a vector 
meson field respectively is also important in reproducing the AA interaction | l8| but in this work we restrict ourselves 
to the most common cr, lo and p mesons. 

The lagrangian density of the baryonic octet model (BOM) reads: 

^BOM = + ^ mesons ~^ ^leptonsi (1^) 

where 

Cb=Y^ ^'b [7p - QvbV^" - gpBt • b'') - {Mb - gsB^)] V'b, 

B 
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with J2b extending over the eight baryons, 

QsB = XsB Qs, QvB = XvB Qv, QpB = XpB Qp 

and XsB, XvB and XpB are equal to 1 for the nucleons and acquire different values in different parametrizations for 
the other baryons, 

i^mesons = ^{9 p,<l)d'^ (j) - mlcj)'^) - ^^Kcj)^ - ^A^^ 

-iB^.•B''- + imX•b^ (16) 

where 12^,^ ~ d^V,y — d„V^, B^^ = 9^b,y — S^b^ — 5p(b^ x b^) and t is the isospin operator. 

In the above lagrangian. neither pions nor kaons are inchidcd because they vanish in the mean field approximation 
which is used in the present worlf and we do not consider the possible contribution of pion and kaon condensates. The 
electromagnetic field contribution also vanishes for homogeneous matter and its effect will not be taken into account 
in the mixed phase. Finally, 

l^leptons = X] {iltid^ - m) tpl- (17) 



In the mean field approximation, the meson equations of motion read: 

B 



6o = ^XpB tsB PB, (20) 

b'^p 

with 

^"" = ^1 (§3(/^+-/^-)' (21) 

1 f M* 
PsB = ^ / P^dp ^ ^={Jb+ + /s-), 



with M% = Mb - 9sB (p, E*{p) = Vp^ + M*^ and /s± = 1/{1 + exp[(£;*(p) tz^b)/^]} , where the effective chemical 
potential is vb = Pb - gvBVo - QpB Ub bo. 

The energy density in the mean field approximation reads: 

-Y^o + Q<f>o + Y<l>o + -Y^o + -Y^o + -Y^o 

+^T.J (^\/p^+W(/^+ +fl-)^ (22) 
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with the leptons distribution functions given after equation H14() . The pressure becomes 



2 ^° 6 24 ■ 2 ■ 24 ■ 2 
+ ^E/^^^(/'+ + /-)- (23) 
The expressions given in this section were obtained for finite temperature, but they can be triviaUy modified for T = 

A. Hadronic matter in beta equilibrium 

The condition of chemical equilibrium is also imposed through the two independent chemical potentials (/i„ and 
/ie) and it implies that: 

Ms" = A^HO = fJ-A = A^n, 

Ms- = A*H- = + /"e, 



For charge neutrality, we must have 



51 toPB + J2 I^P' = 0' (24) 



where qB and stand respectively for the electric charges of baryons and leptons. 

For the hadron phase we have used different parametrizations of the non-linear Walecka model. Two of them have 
been proposed to describe nuclei ground-state properties, NL3 '2^ and TMl |2ll|. and a third parametrization GL 
(in the sequel) was proposed to describe the equation of state of neutron stars p. As shown in the expressions 
given above, the baryonic octet was included. In order to fix the meson-hyperon coupling constants we have used 
several choices discussed in the literature 0, : a) according to 0, we choose the hyperon coupling constants 
constrained by the binding of the A hyperon in nuclear matter, hypernuclear levels and neutron star masses [x„ = 0.7 
and ^ Xp = 0.783) and assume that the couplings to the E and !B are equal to those of the A hyperon; b) we take 
XsB — XyB — XpB = \/2/3 as in 0, H^, 0, 113 ■ This choice is based on quark counting arguments; c) we consider 
XsB — XvB — XpB = 1 known as universal coupling for comparison [2^ . 

In the results we display, whenever Xsb = x^b = Xps, the unique coupling constant will be called xh- We have 
checked that the parametrizations NL3 and TMl become unstable respectively for p/po = 3.4 and 6.5 if xb = ^72/3 
and p/po = 3.5 and 6.7 if xh = 1 owing to the fact that nucleon mass goes to zero. The problem of the effective 
masses turning negative has been discussed in j27l |. and is due to the fact that the model includes baryons with 
different masses and different coupling constants to the a field. This behavior is not exhibited if the hyperons are 
not included, when the effective mass of the nucleons decreases slowly with the density to zero at infinity. However, 
with the inclusion of the hyperons, and if the hyperon couplings are not chosen adequately |27|. the a field grows 
much faster with density. For this reason, we have chosen one of the parametrizations given in |5| (GL) in order to 
obtain the results for the mixed phase discussed in the next section. The chosen parameters are g1/m?g = 11.79 fm^, 
£,2/^2 ^ y -^43 f^2^ gl/ml = 4.41 fm^, n/M = 0.005896 and A = -0.0006426, for which the binding energy is -16.3 
MeV at the saturation density po = 0.153 fm~^, the symmetry coefficient is 32.5 MeV, the compression modulus is 
300 MeV and the effective mass is 0.7Af, higher than in the other two parametrizations. 

In reference [l^ . the authors have imposed a positive value for the effective mass by taking its modulus and 
ignoring the fact that it becomes zero and then negative. This gave rise to discontinuities in their figures for some of 
the parametrizations used. For the TMl parameter set, there are no discontinuities, but in their model they include 
other mesons than the cr, uj and p, and the behavior of the effective mass also depends on the hyperon-meson coupling 
constants used. In ^3] the authors also claim to have used the TMl parametrization in order to study the mixed 
stars, but since no data are given for this parameter set we believe it has only been used for stars with no phase 
transition to the quark phase. 
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IV. MIXED PHASE 

In the mixed phase charge neutrality is not imposed locally but only globally, as mentioned in the introduction 
[slllSjI. This means that the quark and hadronic phases are not neutral separately, but rather, the system will prefer 
to rearrange itself so that 

where is the charge density of the phase i, x is the volume fraction occupied by the quark phase , (1 — x) is the 
volume fraction occupied by the hadron phase and p^. is the electric charge density of leptons. We consider a uniform 
background of leptons in the mixed phase since Coulomb interaction has not been taken into account. According to 
the Gibb's conditions for phase coexistence, the baryon chemical potentials, temperatures and pressures have to be 
identical in both phases, i.e., 

f^HP ~ Pqp, 
Thp — Tqp, 

Php{php,T) = Pqp{pqp,T), 

reflecting the needs of chemical, thermal and mechanical equilibrium, respectively. As a consequence, the energy 
density and total baryon density in the mixed phase read: 

< f X ^"^^ + (1 - X)S"'' + (25) 

and 

<P>-XP'^^ + (l-xK^. (26) 

At this point the EOS in the mixed phase has to be built. In this work we use for the hadronic phase the 
parametrization corresponding io K — 300 MeV and M* /M = 0.7 displayed in the previous section. It predicts a 
mixed phase in a range of densities which depends on the model used for the quark phase and on the hyperon coupling 
constants. The higher the hyperon couplings and the compressibility and the lower the effective mass the lower lies 
the range of densities corresponding to the mixed phase. This has already been discussed in 0,^3- ^'^'^ the quark 
phase we consider both the NJL model and the MIT Bag model with two different bag pressures, as mentioned in 
section II. The choice of the bag pressures was done in such a way that for the lower value, i?^/^ =180 MeV, the onset 
of a quark-hadron mixed phase occurs before the the appearance of hyperons, and for B^/^ = 190 MeV the hyperons 
appear before the quarks. With the present study we have concluded that the existence of a quark phase, and a 
mixed phase is very sensitive to the choice of the hyperon couplings. For a quark phase described by the NJL model 
the mixed/quark phase only exists if the A — coupling is not very weak, namely Xa > 0.65 with the other hyperon 
coupling constants constrained by the binding of the Lambda hyperon in nuclear matter. This can be understood in 
the following way: there will be a quark component only if its effect is to soften the EOS. Large hyperon couplings 
mean that the hadronic EOS becomes harder since at high densities the EOS is dominated by the repulsion described 
by the vector mesons. Therefore a transition to a quark phase may be favorable if it softens the EOS. A weaker Xa 
than a also favors the phase transition for the same reason, i.e., it gives rise to a harder EOS. 

In fig Owe plot the hadronic EOS obtained for several hyperon couplings and the NJL EOS as a function of the 
neutron chemical potential For the hyperon couplings we have considered besides the couplings already introduced 
in the present paper, two other sets, both constrained by the A hyperon binding to nuclear matter with Xa = 0.6 
and 0.8. For the first value of Xg- we do not get a phase transition to quark matter with the NJL model. The value 
x„ = 0.8 was chosen in ref. 14], but according to 5] could be a bit too high if we consider hypernuclear levels 
which give an uncertain upper bound of x„ = 0.72. The EOS with xh=^- is very similar to the one obtained with 
Xa = 0.7, Xi^ = Xp = 0.783. We have assumed that the couplings to the S and S hyperons are equal to those of the A 
hyperon. However, we have also tested that if we had taken for the S couplings the two different sets proposed in [2^ . 
which give satisfactory fits to the S" atom (xaT. — 0.77, x^iT. — 1-, XpY, = 0.67 or XaY, — 0.54, Xuiy: — 0.67, Xps — 0.67) 
the densities corresponding to the onset of the mixed phase and the pure quark matter phase would not alter. 

According to for a newborn neutron star the entropy per baryon across the star is approximately constant, of 
the order 1 to 2. This corresponds to temperatures which can go up to 30 MeV In fig 21 we plot the entropy 

for different temperatures obtained with the NJL model. We conclude that, except, for a small effect in the mixed 
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phase, in the interior of the star {p > 2po) the entropy does not change with temperature. Exactly the same trend is 
observed with the Bag model. Therefore, we use the EOS obtained for different parametrizations at fixed temperatures 
to determine the properties of proto- neutron stars, and we are confident that there will not be big differences for not 
using fixed entropies. 

The rise of the entropy in the mixed phase corresponds to a decrease in temperature in the interior of the star, if a 
fixed entropy calculation were performed. This effect has also been discussed in ref. |Q| and it is due to the opening 
of new degrees of freedom in the system. Similar results were obtained with the Bag model. 

In order to implement both phases in a common code, we have imposed the Gibb's conditions and rewritten the 
quark chemical potentials of equations (|I3|I in terms of the electron and neutron chemical potentials, i.e., 

^J■u = {^J■7l - 2Aie)/3 Pd= tJ'S = iPn + A*e)/3. 

Examples of the particle populations obtained from the construction of the full EOS are shown in figures El and 
Notice that if the Bag model is used, the quarks appear at relatively low densities, sometimes lower than the densities 
for which some of the hyperons turn up, depending on the value of Bag. This behavior was also observed in P, |^. 
If the NJL model is used the appearance of the quarks depends very much on the hyperon coupling constants. Two 
different behaviors were obtained: a) the hyperons appear before the quarks and the mixed phase occurs at much 
higher densities, (~ 5po); b) the quarks appear at T = at lower densities than the hyperons if the effective bag 
pressure at the onset of the mixed phase is still well below the plateau around 162 MeV. In this second case the phase 
transition occurs at ~ 2po and the effective bag pressure increases to its plateau value for densities inside the mixed 
phase. This gives rise to several effects such as a plateau on the quark content after the initial increase as can be 
observed in figure[7|at T = 0. For T > 10 MeV this last effect disappears because the rise of the effective bag pressure 
occurs at higher densities. With this set of hyperon couplings there is a transition with temperature: for T = MeV 
the quarks appear before the hyperons, at T = 5 MeV both appear at a similar density; for higher temperatures the 
hyperons appear at lower densities. One main difference between the Bag model and the NJL model predictions is 
the relative fractions of u, d, and s quarks: in the Bag model all the three types of quarks appear at a similar density, 
then the fraction of d and s quarks increases faster and only close to the onset of the pure quark matter the u fraction 
becomes larger than the s quark fraction. At these densities, however, the u, d, s quark fractions are very close to 
each other, of the order ~ 0.33. With the NJL model the d quark appears first, at a higher density the u quark 
appears and still higher the s quark turns up. The d{s) quark fraction is always the largest (lowest). Only at very 
high densities, well inside the pure quark matter phase, do all three quark fractions get close to 0.33. These behavior 
is due to the fact that in the NJL model, only at quite large densities the mass of the s quark gets close to its current 
mass value. 

For all EOS obtained the higher the temperature, the lower the densities for which the hyperons appear. The onset 
of the mixed phase occurs in most cases studied at higher densities for a higher temperature, the exception being the 
the NJL model for xh = ■\/2/3. This can be confirmed in Fig. [3 where the full EOS, obtained respectively with the 
Bag and the NJL models, for xh — ■\/2/3, are given for different temperatures. For this choice of hyperon couplings 
the Bag model predicts a mixed phase at lower densities (energies) . For the other two choices of hyperon couplings 
studied with the NJL model, the mixed phase occurs before the effective bag pressure reaches the 163 MeV plateau, 
and this dictates the difference. For Xs = 0.7 and xh = 1-0 the onset of the mixed phase is determined by the low 
density behavior of the NJL EOS. As discussed in j23|, the EOS obtained within this parametrization of NJL model 
contains mechanical unstable regions which give the possibility of a first-order phase transition of the liquid-gas type. 
Depending on the set of hyperon couplings chosen, it can happen that the mixed phase appears precisely at the 
density above which the NJL is stable against density and charge fluctuations. 

The presence of strangeness in the core and crust of neutron and proto-neutron stars can have important conse- 
quences in understanding some of their properties [2^180. .31. 32]. We have calculated the strangeness content of the 
different EOS as a function of density. In figures |H1 and we plot the strangeness fraction defined as 

rs=xr?'' + {l-x)rr (27) 

with 

OP ^ Pf. HP ^ J2b lis \Pb 
3p' ' 3p 

where qf is the strange charge of baryon B. When the quark phase is described by the Bag model the strangeness 
fraction rises steadly and at the onset of the quark phase it has almost reached 1/3 of the baryonic matter. This 
behavior is independent of the hyperon-meson coupling constants used in this work. The NJL model predicts a 
different behavior: in the mixed phase the strangeness fraction decreases and increases again for pure quark matter. 
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This behavior is due to fact that for the densities at which the mixed phase occurs the mass of the strange quark is 
still very high. The overall effect of temperature is to increase the strangeness fraction, except for the mixed phase 
with the NJL model where the strangeness fraction decreases more strongly for higher temperatures. 

V. MIXED STAR PROPERTIES 

In the present section we investigate the properties of stellar objects formed by matter described by the EOS studied 
in the previous sections. Mixed proto-neutron and neutron star profiles can be obtained from all the EOS studied by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equations _33J, resulting from a simplification of Einstein's general 
relativity equations for spherically symmetric and static stars. For a certain EOS the TOV equations for pressure and 
mass are integrated from the origin for a set of arbitrary choices for the the central density, in such a way that they 
define a one-parameter family of stars. The neutron star mass is an increasing function of its central density until 
it reaches a maximum, where matter can collapse to a black hole. It means that the values of the maximum mass 
of a star is very important in the determination of the mass of possible black holes. The Oppenheimer-Volkoff limit 
gives the maximum value of the neutron star mass. There are theoretical previsions that the maximum mass varies 
from 0.72 solar mass (Mq) obtained for pure neutron matter stars up to 3.26 M© 0|. Other authors based on 
recent calculations, claim that they can be in the interval 1.6 Mq and 2.1 Mq. According to the observed values 
lie between 1.2 Mq and 1.8 M©. 

In tables 1 and 2 we show the values obtained for the maximum mass of a neutron star as function of the central 
density for some of the EOS studied in this work and for different temperatures. In table 1 the GL force and the NJL 
model were used to derive the full EOS, while in table 2 the EOS was obtained from the GL force and the MIT Bag 
model. Several conclusions can be drawn. For all EOS studied, the central energy density Eq falls inside the mixed 
phase. In some cases, at finite temperature, it can even occur in the quark phase. This happens in the calculation 
with NJL if the universal coupling is chosen and in the calculation with the Bag model with B-^^'^ = 180 MeV. Similar 
results were obtained in [31 ■ However, the present results are different from the ones discussed in [l^, where the 
parametrization of the EOS used allowed for the appearance of the mixed phase only in a quite narrow mass range. 
The maximum masses of the stars do not show any regular behavior with temperature, contrary to the steady increase 
of the maximum allowed mass obtained with the GL model alone as shown in table 3. It is worth emphasizing that 
the EOS for T=0 and T=10 MeV are very close to each other (at least for GL plus Bag models), what could have 
led to very small differences in the solutions of the TOV equations. The trend seems to be a small decrease at lower 
temperatures followed by a small increase. Only in one case we have obtained a continuous increase of the mass. A 
discussion of the possibility of a back hole formation after the deleptonization and cooling of the proto-neutron star 
requires a more careful study. In principle, this could happen if the maximum allowed mass at T = MeV is lower 
than the maximum allowed mass at finite T. 

Comparing tables 1 and 2 we conclude that the Bag model allows for smaller maximum masses, of the order 
~ 1.6 Mq, than the NJL model, ^ 1.9 Mq. Several quite high mass limits have been recently determined [ssl l3^ 
which could rule out soft EOS with phase transition to deconfined matter. In particular, the determination of the 
pulsar Vela X-1 gives the mass Mx — 1.86 ± 0.16 Mq. However, these high mass limits do not exclude a phase 
transition to quark matter within the NJL model. It is also seen that the maximum mass does not depend much on 
the hypcron couplings although in the NJL the same is not true for the onset of the mixed phase. For comparison we 
show in table 3 the properties of stars obtained with the GL model (transition to the quark phase not included). While 
the maximum masses obtained with Bag model are ^ 15 % smaller than with the hadron EOS, the ones obtained 
within the NJL model are only 3% smaller. 

We will not discuss the radius of the maximum mass star because it is sensitive to the low density EOS and we did 
not describe properly this range of energy densities. 

VI. CONCLUSIONS 

In the present paper we have studied the EOS for proto-neutron stars using both the Bag model and the NJL model 
for describing the quark phase and a relativistic mean-field description in which baryons interact via the exchange 
of a—, oj—, p— mesons for the hadron phase. Parametrizations TMl and NL3 fitted to the ground-state properties 
of stable and unstable nuclei proved to be inadequate because, due to the inclusion of hyperons, the nucleon mass 
becomes negative at relatively low densities. We have considered a parametrization which describes the properties of 
saturating nuclear matter proposed in ^ . For the hyperon couplings we have used three choices and verified that for 
the NJL model the onset of the mixed phase is sensitive to the hypcron couplings. We have verified that the choice 
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of other sets of hyperon couplings, namely couplings to the E and 5 different from the A hyperon, did not have any 
effect on the onset of the mixed phase. 

If the quark phase is described by the Bag model the strangeness content of the EOS rises steadly and at the onset 
of the quark phase it has almost reached 1/3 of the baryonic matter. However, the NJL model predicts a different 
behavior: in the mixed phase the strangeness fraction decreases and increases again for pure quark matter. The overall 
effect of temperature is to increase the strangeness fraction, except for the mixed phase with the NJL model where 
the strangeness fraction decreases more strongly for higher temperatures. The strangeness content will influence the 
cooling rates of the star 31, 32]. 

The consequences of using a model with chiral symmetry have been discussed. The onset of the mixed phase is 
the result of a delicate equilibrium between the hyperon couplings and the rise of the effective bag pressure due to 
partial restoration of the chiral symmetry. The value of the effective bag pressure also dictates whether the onset of 
hyperons occurs at larger or lower densities than the onset of the mixed phase. Temperature will alter the order if at 
T = MeV the onset of the mixed phase occurs at lower densities. A large s quark mass until very high density gives 
rise to a decrease of the strangeness content of the EOS in the mixed phase. In the quark phase the fraction of the s 
quark becomes similar to the d and u quark fractions at densities greater then 10 po- 

The appearance of the mixed phase in /3-equilibrium matter at a fixed temperature gives rise to an increase of the 
entropy per baryon due to the opening of other allowed channels. This will affect neutrino-matter interaction rates 

m 

With the EOS studied in the present work, the mixed stars obtained from the integration of the TOV equations 
contain a core constituted by a mixed phase. Only in very special cases, a small core of pure quark matter exists. It was 
shown that the maximum masses of mixed stars obtained with the Bag model are of the order of '--^ 1.6 Mq , smaller 
than the maximum masses obtained within the NJL model, ~ 1.9Mq. Even with a transition to a deconfinement 
phase the masses predicted can be quite high. Maximum masses for neutron stars obtained in the present work are 
larger than the ones obtained in P|. The effect of temperature in the maximum masses allowed is not strong. The 
central energy densities, though, decrease with temperature. 

We are aware of the importance of neutrino trapping in the description of proto-neutron star properties and this 
feature will be incorporated in a forthcoming work. The effects arising from the inclusion of pion and kaon condensates 
seem not to be very large and also very dependent of the choice of the hyperon to meson coupling constants_ 
which are not known quantities. For these reasons they were not incorporated in this study. Some authors 
have claimed that other meson fields, namely the scalar meson field /o(975) and the vector meson field 0(1020), 
should also be incorporated in order to reproduce the observed strongly attractive AA interaction. We believe that 
the inclusion of these meson fields will not have a great influence on the properties of the stars under consideration. 
It is well known that non-rotating neutron stars are practically undetectable. The extension of this work in order to 
include the rotation mechanism discussed in l35| is under consideration. 

Moreover, the quark-meson coupling model (QMC) [s^, which describes nuclear matter as a system of non- 
overlapping MIT Bags which interact through the effective scalar and vector mean fields, is also a good option 
for the description of the hadron phase. The quark phase appears naturally as the breakdown of the confining bags 
takes place. This investigation is also under way. In the present work we have only considered unpaired quark 
matter. However, recent studies seem to show that the quark phase, if present, would be in a color superconductor 
phase [39I l40| . The possible existence of a color superconductor quark phase in compact stars has been studied in 
references |4ll | with the Bag model and using the NJL model. In general, it was shown that the inclusion of a 
color superconductor phase would soften the EOS at lower pressures giving rise to the onset of quarks, either through 
a sharp transition or a mixed phase, at lower densities. With the NJL model and a sharp transition assumption no 
stable stars with a quark core were found. Properties of compact stars obtained within the Bag model depend on the 
parameter set used and the sharp/mixed transition considered, but it was shown that color superconductivity could 
boost the mass of the star. Further studies on the possible existence and consequences of a color superconductor 
phase are still needed. 
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Table 1 - Mixed star properties for the EOS obtained with the GL force and the NJL model. 



T (MeV) Mmax/Mp ep (fm-^) Smin (fm"^) Sma. (fm-^) 



XH = vVs 





1.84 


6.29 


4.60 


7.25 




10 


1.83 


6.34 


4.58 


7.14 




20 


1.84 


6.26 


4.50 


6.66 




30 


1.85 


5.85 


4.23 


5.84 


Xs = 0.7 





1.90 


5.01 


1.30 


5.21 


x^ = 0.783 


10 


1.89 


4.98 


1.31 


5.13 


X(^ = Xp 


20 


1.89 


4.81 


1.37 


4.82 




30 


1.90 


4.34 


1.61 


4.44 


XH = 1. 





1.92 


5.08 


1.30 


5.22 




10. 


1.89 


5.13 


1.54 


5.10 




20 


1.88 


4.92 


2.20 


4.83 




30 


1.87 


4.58 


2.18 


4.46 



Table 2 - Mixed star properties for the EOS obtained with the GL force and the MIT Bag model. 



T (MeV) Mmax/Mp eg (frn'^) (fm-^) Smax (fm-^) 



BagV*=180 MeV 





1.40 


7.38 


1.17 


4.62 


XH = 72/3 


10 


1.39 


7.05 


1.17 


4.58 




20 


1.42 


4.18 


1.20 


4.47 


Bagi/4=190 MeV 





1.64 


4.58 


1.81 


6.06 


xh = 1/2/3 


10 


1.59 


4.86 


1.85 


6.03 




20 


1.67 


4.18 


1.90 


5.93 




30 


1.76 


3.18 


1.96 


5.78 


Xs = 0.7 





1.63 


4.43 


1.53 


6.0 


x^ = 0.783 


10 


1.63 


4.41 


1.57 


5.95 




20 


1.63 


4.10 


1.62 


5.86 




30 


1.64 


3.82 


1.67 


5.69 


Xh = 1 





1.64 


4.49 


1.63 


6.01 




10 


1.65 


4.40 


1.66 


5.98 




20 


1.67 


4.12 


1.71 


5.87 




30 


1.72 


3.61 


1.76 


5.70 



Table 3 - Compact star properties for the EOS obtained with the GL force 
T (MeV) Mma^/Mp eg (fm"^ 



XH = ^/273 






1.93 


6.36 


10 


1.93 


6.34 


20 


1.95 


6.20 


30 


1.97 


5.85 




FIG. 1: NJL model in beta equilibrium: the bag effective pressure at different temperatures 
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FIG. 2: The fraction of quarks u,d,s in beta equilibrium at T=0, 20 and 50 MeV (from very thick to thin) 




FIG. 3: The EOS for the NJL model and for GL model with different choices of the hyperon couplings) 
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FIG. 4: Entropy per baryon for different temperatures. The quark phase is described by the NJL model. 
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FIG. 5: EOS obtained with the GL force plus a) Bag model for Bag^''*=190 MeV; b) NJL model. In both cases xh = 
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FIG. 6: Particle fractions Yi = Pi/p, for i = baryons, leptons and quarks, obtained with the GL force plus the Bag model 
Bag^/*=190 MeV and a) T=0 b) T=20 MeV. In both cases xh = x/273. 
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FIG. 8: Strangeness fraction for the EOS with the Bag model for the quark phase with B^^^ = 190 MeV and xh = sjlj?). 




